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Beyond adiabatic elimination: Effective 
Hamiltonians and singular perturbation 


Mikel Sanz, Enrique Solano, and Inigo L. Egusquiza 


Abstract Adiabatic elimination is a standard tool in quantum optics, that produces 
an effective Hamiltonian for a relevant subspace of states, incorporating effects of its 
coupling to states with much higher unperturbed energy. It shares with techniques 
from other fields the emphasis on the existence of widely separated scales. Given 
this fact, the question arises whether it is feasible to improve on the adiabatic ap¬ 
proximation, similarly to some of those other approaches. A number of authors have 
addressed the issue from the quantum optics/atomic physics perspective, and have 
run into the issue of non-hermiticity of the effective Hamiltonian improved beyond 
the adiabatic approximation. Even though non-hermitian Hamiltonians are interest¬ 
ing in their own right, this poses conceptual and practical problems. 

Here, we first briefly survey methods present in the physics literature. Next we 
rewrite the problems addressed by the adiabatic elimination technique to make ap¬ 
parent the fact that they are singular perturbation problems from the point of view of 
dynamical systems. We apply the invariant manifold method for singular perturba¬ 
tion problems to this case, and show that this method produces the equation named 
after Bloch in nuclear physics. Given the wide separation of scales, it becomes intu¬ 
itive that the Bloch equation admits iterative/perturbative solutions. We show, using 
a fixed point theorem, that indeed the iteration converges to a perturbative solution 
that produces in turn an exact Hamiltonian for the relevant subspace. We propose 
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thus several sequences of effective Hamiltonians, starting with the adiabatic elimi¬ 
nation and improving on it. We show the origin of the non-hermiticity, and that it is 
inessential given the isospectrality of the effective non-hermitian operator and a cor¬ 
responding effective hermitian operator, which we build. We propose an application 
of the introduced techniques to periodic Hamiltonians. 


1 Adiabatic Elimination 

Situations in which there is a wide separation among energy or time scales present 
in a physical system pervade all of Physics. It is therefore of recurrent interest to de¬ 
velop techniques to obtain approximate, effective descriptions of the low energy or 
slow sector of the system, since they are the most likely accessible to experimenta¬ 
tion and control. In general, a naive perturbation expansion will not provide us with 
the required effective description; depending on the scheme and approach, the obsta¬ 
cle to do so will turn up as secular terms in a time evolution, or as zero denominators 
in state expansions. Therefore, the approximation schemes valid for these situations 
will not be directly perturbative; they will be asymptotic, or resummation-based, or 
cumulant, or combinations thereof. Often, we have heuristic arguments for the con¬ 
struction of an effective description of the system, which do not naturally lead to 
improvements, at least not systematically. 

A very common and useful approximation in quantum optics, normally based on 
a heuristic argument, is the so-called adiabatic elimination technique EL®®. There 
are several ways of introducing this approximation. For this first presentation, con¬ 
sider the Schrodinger equation as a dynamical system, evolving with a Hamiltonian 
H as id t \f/ = Hy/ (we set here and henceforth Ti = 1). Let state y/ be partitioned into 
a = Py/ and y = Qyr, with P and Q = 1 — P projectors, in such a manner that the 
eigenvalues of PHP are widely separated from those of QHQ. Let it be the case in 
which the coupling between the P and Q subspaces is very small when compared to 
the eigenvalues of QHQ. To be more specific, let T be a characteristic scale, for in¬ 
stance, the norm of the restricted inverse {QHQ)~ l . Then, the Schrodinger equation 
may be rewritten as 


id, a = PHP a + PHQy , 
i d,y= QHPa + QHQy , 


(la) 

(lb) 


where the last line of the system is to be multiplied by T. Heuristic ally, we are ask¬ 
ing that T QHP <C 1, while T QHQ = (9(1). This can be achieved if y is small, and 
if we neglect its time evolution. In other words, we slave y to a in the approxima¬ 
tion QHQy= —QHPa , or, formally, y= —{QHQ)~ l QHPa. By substituting this 
approximation for y in the first component of the system, we obtain the effective 
evolution for the slow sector a as 



(2) 
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By eliminating the fast component 7 , we obtain an effective Hamiltonian for the 
slow sector 

//adiabatic = PHP - PHQ-^—QHP . (3) 

This adiabatic elimination process has proven to be extremely useful in isolating 
the effective slow evolution for a number of systems . However, the aforementioned 
heuristic presentation does not lend itself readily to a systematic improvement. Thus, 
the question of how to transform this approximation into a controlled expansion is 
recurrently posed in quantum optics literature 00 . In essence, both papers present 
an expansion of an integral kernel, either in the energy domain 0 or the time do¬ 
main Q. 

As a matter of fact, this problem of finding a systematic expansion when there 
is a wide divergence in energy scales has also appeared in many other contexts: 
nuclear physics, condensed matter, and atomic and molecular physics, among oth¬ 
ers. In each case, there are expansion techniques, under different guises and names: 
Bloch equation J 6 | 0, Schrieffer-Wolff expansion (S), or Born-Oppenheimer ap¬ 
proximation 0 . In this paper, we want to merge these techniques by presenting 
the system of equations as a singular perturbation problem, and by using the invari¬ 
ant manifold scheme for resummation of secular terms. For this purpose, we shall 
derive systematic expansions which improve on the adiabatic elimination technique 
and furnish us with an effective Hamiltonian for the low energy sector. Furthermore, 
we show the equivalence of this expansion and the Schrieffer-Wolff method. We ap¬ 
ply this method to some simple examples in quantum optics, and suggest its future 
applicability for time-dependent periodic Hamiltonians. 


2 A singular perturbation problem and Bloch’s equations 

In Physics, there are often problems depending on a small parameter which can be 
solved by applying perturbation theory. The singular perturbation theory approach 
has been developed for those cases for which no uniform regular expansion is pos¬ 
sible, the so-called singular problems. Heuristically, one can identify a problem as 
being singular when it is qualitatively different at the zero order of the expansion, 
and immediately out of it. For example, the algebraic problem ex 2 +x — 1 = 0 is 
singular for the small parameter e, since it is second order if £ ^ 0 , but first order if 
£ = 0. Similarly, singular perturbation problems in systems of ordinary differential 
equations m appear frequently by reducing the order of the differential equation 
when the small parameter is set to be 0 . 

Let us now consider the A system depicted in Fig. [T] which is ubiquitous in 
atomic physics. This is governed by the system of ordinary differential equations 
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Fig. 1 The A system is a paradigmatic model in atomic physics and quantum optics, where a 
three-level atomic system may be coupled to classical or quantum fields. 


8 £ 2 * 

ia = -2 a + ^ 7 ’ 


i/3 




(4) 


■ . n , 

i7= ^-a + ^-/3+4y. 


We are interested in the regime in which 72, |, with 5 real. Under the change 

of variable t 8t and the definitions Q, = Q,/8 and e = 8/A, the system trans¬ 
forms into 


1 £ 2 * 

m = --a + — r , 




O-: 


-7, 


(5) 


iey= y + -(Q a a + Q b li) . 


This is heuristically identifiable as a singular problem, since the last differential 
equation changes its character into an algebraic one as £ —> 0. Indeed, if one were to 
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attempt a naive perturbation expansion in the small parameter e, one would obtain 
secular terms already at the first order, rendering the expansion invalid already for 
times t of the order of A /| A (in terms of the original variables). 

The issue at hand is therefore how to identify a time-uniform scheme which pro¬ 
vides us with approximations to the slowly varying variables. Among the numerous 
existing proposals, we shall concentrate here on the invariant manifold method iflOl . 
This consists in constructing, among the submanifolds invariant under the flow, 
those that are perturbative in the small parameter. In the case of interest, as the 
system is linear, it is pertinent to only examine linear subspaces. 

Let us write system G in the matrix form 


id t 





(6) 


For clarity, let us notice that if we were to write the A system in this more compact 
notation, we would be using the isomorphism C 3 = C 2 © C. The three component 
vector (a,j3,y) r would be written (again for the A system) as a two component 
object (a,y) T , where the first component is in turn a vector in C 2 . It follows that, 
in this case, ft) £ Mi(C) (in fact, hermitian), A £ M i(C) and hermitian, so in fact 
A £ K. for the A system case. Finally, Q is a row two component vector for the A 
system. 

For the general case, consider that the full Hamiltonian of system G acts on 
the Hilbert space M. Then the projectors P and Q give the Hilbert subspaces PM 
and QM respectively, with PM' © QM' = M' . We have introduce new notation 
in (|6ji for the operators PiHPj, where P, stand for either P or Q. Namely, PHP — > 
co £S$ (PM), QHQ —> A e 3S{ QM ’), while QHP -)• Q £ B§{PM’,QM’). Now, a 
stands for an element of PM, while y £ QM 1 . 

As we are looking for linearly invariant subspaces, let us define them by the 
relation y = Ba. Notice that the embedding operator B belongs to 34{PM\ QM r ), 
as Q. The invariance condition or embedding equation then reads 

Q +AB = B(0+BQ^B. ( 7 ) 


By itself, this equation, known in the literature as Bloch’s equation nn, is not ad¬ 
vantageous with respect to the direct analysis of the Hamiltonian. However, in the 
case of interest to us, the spectrum of A is energetically very separated from the 
spectrum of ft), and we shall carry out a perturbative expansion or an iterative pro¬ 
cedure to determine B. Once B has been obtained, the evolution in the subspace is 
determined by i d t a = (ft) + £> B) a. This also provides us with approximate evo¬ 
lutions in the subspace, given an approximate solution to Bloch’s equation. Notice 
the existence of an, in general, non-hermitian linear operator h e g- = ft) + Q ' B, which 
plays the role of an effective Hamiltonian. 
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3 Expansion beyond adiabatic elimination 


Let us assume that A has a bounded inverse, and that, for some definition of the op¬ 
erator norm, e = ||A~ 1 ||||a)|j <C 1 and e' = ||A _1 || ||f2|| -C 1. We define the nonlinear 
transformation of operators 

T(A) = -A- l Q+A~ l Aco + A- l AQ f A. (8) 

The invariance condition (j7]i may be now written as a fixed point equation, namely 
B = T ( B). We shall now prove that the nonlinear transformation T has a fixed point 
which is the required solution to the problem at hand. Observe that, due to the defi¬ 
nition of T and the properties of operator norms, we have 

II^WII ^ ll^ _1 || ll^ll + ll^^ 1 |l ll®ll ll^ll + ll^ _I || ll^ll ll^ll 2 

< e'(l + ||A|| 2 )+e||A||. (9) 

From this observation we obtain the following central proposition: 

Theorem 1. Assume e, e' > 0 and e' < (1 — e)/2. Let us define 


- 1 - 


( 10 ) 


Then, the fixed point equation T{A) = A has at least one solution A* such that 


||A*|| < r(e,e'). (11) 

Proof. By direct analysis of the function g(x) = e'(l +v 2 ) + £x, and using inequality 
([9]) we conclude that, for every A such that 

l|AII < r(e,e') (12) 


it holds that 

||r(A)||<r(e, e '). 03) 

As T maps a bounded closed convex set of the corresponding Banach space of 
operators into itself, there exists at least one fixed point A* = T(Af) in the set, by 
Schauder’s fixed point theorem. 

Furthermore, choosing the operator — as the initial point of the iteration, 

one can readily see that 

II—A- 1 f2||<||A- 1 ||||f2||= e '<r(£, e ') (14) 


under the stated conditions. Thus, we know that there is at least one fixed point in 
its vicinity. 

To sum up, under those conditions, we are assured of the existence of a solution 
of Bloch’s equation which is small in the sense that ||B|j < r(£,£'). It is therefore 
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natural to attempt either an iterative process or a perturbative expansion to compute 
approximations for that solution. 

More concretely, we define 


fi (0) = -A~‘Q, 


_ J 



(15) 


The sequence of operators B lk ) is thus defined by iteration, and lies in the region 
of applicability of the proposition. Although we have not proven the convergence 
of this sequence, the mapping of problem ([7]i to the Schrieffer-Wolff expansion 
presented in Sec. fallows us to conclude the uniqueness of the fixed point in some 
circumstances, given the results presented in fT2i . 

Alternatively, we can define a perturbative expansion as 

B=i B ( k) < 16 ) 

k—l 

and insert it into the embedding or Bloch equation to obtain the recursive equation 

k— 1 

B( k +i) =A- l B [k] co + A~ l (17) 

/= 1 

with the initial conditions Bm = -4 1 Q and Bro) = —4 2 Qco. Notice that, for¬ 
mally, 

fi(fc) -L fi (0=°(4" ( * +2) ) . 

i=t v 7 

which relates the perturbative and the recurrence results to a given order. Explicit 
formulae relating the iteration approximations B^> and the perturbative contribu¬ 
tions Sm are hardly illuminating beyond this relation to a given order. 


4 The effective Hamiltonian 

As pointed out above, the evolution in the subspace fixed by B is determined by 
ft) I Q B. In the first order, this provides us with a hermitian operator 

h w =/;(!) = co-Q J( A-'n, (18) 

either by iteration or perturbatively. This is, in fact, what we would obtain from 
direct adiabatic elimination, namely, by setting <9,y to zero, solving y as —4 
and replacing this for a in the differential equation. Furthermore, this operator is 
hermitian by construction. 
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Hermiticity, however, is not maintained in higher recurrence or perturbation or¬ 
ders ED. This has been a sticking point in the literature, and the source of some 
confusion. It is clear that, by construction, the linear generator of evolution for the 
a part of the full t ff = (a, y) 7 does not need to be hermitian, even if the total Hamil¬ 
tonian is hermitian. In any case, its spectrum must be real. In order to prove it, let 
us consider an eigenvector of the total Hamiltonian belonging to the linear subspace 
determined by y = Ba. Then, the restriction to its a part will be an eigenvector of 
the effective Hamiltonian with the same eigenvalue as that for the total Hamilto¬ 
nian. Conversely, let us assume that a* is an eigenstate of the effective Hamiltonian. 
Then, y/* = (ct*,Ba*) r is also an eigenstate of the total Hamiltonian with the same 
eigenvalue, and hence the latter is real. 

In the finite-dimensional case of interest, it follows that the full effective Hamil¬ 
tonian must be similar to a hermitian Hamiltonian. Let us construct this similarity 
transformation assuming that we already have a solution B for Eq. 0 at hand. As a 
first step, notice that the total conserved norm (ty, yr) can be expressed as 

{¥,¥) = (a,a) + (y,y) 

= (a, (l +B T B) a). 

One is immediately led to examine hi = (l +B'B) (w + 'fi). By using Eq. (j7|i 

and its conjugate, one concludes that h\ is indeed hermitian, if B is a solution to 
Bloch’s equation. Its spectrum, however, is not the one corresponding to the time 
evolution in the low energy sector PJff. Let us define 

S B = Vl+B^B, (19) 

which is always possible, since 1 +B'B is positive. Then, for any constant unitary 
V acting on PJt?, we obtain a hermitian Hamiltonian 

h v = VS b ((o + &B)S b 1 V^ . (20) 

Hermiticity is easily proven by noticing that 

h v = VS B (co + Q Jf B)S B l V t 

= vs B l h x s B l v\ 

The unitary V can be subsumed in the choice of the square root defining Sb in ( p~9| >. 
In what follows, the omission of the subscript y is associated to the assumption of a 
choice for Sb- 

In point of fact, an exact solution for (|7]i is as hard to come by as an exact di- 
agonalisation of the initial Hamiltonian, so we have to use approximate methods. 
However, if we truncated B to some approximation, hi (and thus h v ) would no 
longer be automatically hermitian. This can be mended by introducing in 

h = (l+B t B)“ 1/2 (© + X2 t B + B t f2+B t AB) (l+fi f fi)“ 1/2 


( 21 ) 
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an approximate solution B a to Bloch’s equation. This expression would be equiva¬ 
lent (up to unitary transformations V ) to hy, if B a were an exact solution of Bloch’s 
equation ([7]); on the other hand, it is explicitly hermitian for any B a , which provides 
us with a hermitian approximate Hamiltonian encoding the effective evolution. In 
the case of a perturbative expansion, one can compute m to second order, for 
example, giving 

M 2) = * (&A- 2 na)+(on f A- 2 n). (22) 


5 Relation to Schrieffer-Wolff expansion 


In 1966, J. R. Schrieffer and R A. Wolff introduced a perturbatively-built canon¬ 
ical transformation, with an antihermitian generator S, in order to eliminate small 
denominators in the perturbative expansion of the Anderson Hamiltonian QjQ. How¬ 
ever, this technique is nowadays a staple in condensed matter physics. 

Let us define an antihermitian operator S by 


tanh(S’) = B Q ^ , 


(23) 


where S is called Schrieffer-Wolff’s operator (see, for instance, Ifl2l ). Then, Bloch’s 
equation 0 and its conjugate can be written together in the form 


0 

Q 0 


— tanh(S) 


0 £2 T 
Q 0 


tanh(S) 


(o a ) ’ tanh(5) 


Under the assumption that cosh (S'), with antihermitian and block anti-diagonal S), 
exists and is invertible, this equation is equivalent to exp(5j//exp(— S) being block- 
diagonal. 

Some further formal manipulations give us the additional identification 


e s = [1 — tanh(S)] [l — tanh 2 (S)] 

_ (\ /( l - t - fi 1 '#)' 172 0 

~\ B 1 y V 0 (1 +BB^y l/2 


(24) 


This expression, together with the equivalence between Bloch’s equation and that 
exp(5’)//exp(— S) is block diagonal, give us the equivalence between the two meth¬ 
ods. The hermitian effective Hamiltonian hy is, in fact, the block obtained from the 
Schrieffer-Wolff method corresponding to the low energy sector, and an approxi¬ 
mate S gives us an approximate B, and viceversa. 

In the Schrieffer-Wolff method, the approximations are built by expanding S in 
powers of the interaction, and imposing that, to that order, the transformed Hamilto¬ 
nian must be block-diagonal. To establish the connection between both methods, let 
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( 0 — 

us denote .S' = I ^ ' ). Then, by using our notation, the first order Schrieffer- 

Wolff condition reads 


yb(o-Ayb = n. 


(25) 


Notice that this is a Sylvester equation, which will have a unique solution if the 
spectra of ft) and A are disjoint. It should be pointed out that, in the context we are 
interested in, it might be necessary to compute approximate solutions for Eq. 
Alternatively, if a solution for (25 i can be computed, it will be a resummation of our 
perturbative or iterative expansion. 

To the same expansion order, the effective Hamiltonian for the low energy sector 
reads 

h ( ^ = o)+ ] -^n+n f 


■ 


(26) 


6 The A-system example 


The A-system in quantum optics described in the system of Eqs. 0 is amenable 
to both the approach based on Bloch’s equation and Schrieffer-Wolff’s approach if 
A 8, Direct application of expressions (22 1 and (261 yields 


hsw 2°' Al-{8/2A) 2 


+ 


4 A 2 1 — (8/2A) 2 


o z &Q -t-iiT.Ocr') , 


/z (2 ) =- o' — 

2 A 

C 

+ \ (o : n f Q +&&o z ) . 

4A 2 v ' 


(27) 


(28) 


In this simple example one can already see some features of the Schrieffer-Wolff 
expansion as compared to the Bloch expansion; namely that the coefficients are not 
purely perturbative, but involve a resummation of perturbative terms. Additionally, 
this resummation presents with a pole that is not seen in the first perturbation terms. 
One should not expect the location of the pole to this order to be exact, and, in fact, 
working out specific exactly solvable examples (such as Q. —> (1,0)) it is easy to see 
that it is not located at 8 = 2A 

Alternatively, the iteration expansion readily lends itself to numerical implemen¬ 
tations. In Fig. [2] and [3] we show the presence of a secular shift for the adiabatic 
elimination approximation, the lack of normalisation with the non-hermitian Hamil¬ 
tonian to some orders, and that the numerical solution of Bloch’s equation matches 
perfectly to the real one already at fourth order of iteration, when ignoring large 
frequency oscillations. 
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Fig. 2 Evolution of the population of the ground and excited states, with initial state (l,0,0) r , 
under a) the exact Hamiltonian (continuous black line), b) zeroth order effective Hamiltonian 
(dashed blue line) and c) fourth iteration of T (dotted red line). The parameters are S = —0.01754, 
Q a = 0.44, Qb = 0.34, for direct comparison with Ref. (5j. 


In particular, in Fig. [2} we show the exact evolution of the populations of the 
lowest lying states in the A system. In the context of atomic physics, “populations” 
translates into the square modulus of the coefficients of the state in an orthogonal 
basis. 

Thus, for a state l// = (a,p,y) T we are depicting |a| 2 and |/3| 2 as a function of 
normalised time At, and these two quantities are compared with the evolution of the 
same populations with the first order of Bloch’s approximation, i.e. adiabatic elim¬ 
ination approximation, and with the fourth iteration of the recurrence. One should 
observe the secular shift in the adiabatic elimination: the maxima of the relevant 
population under the evolution dictated by the adiabatic elimination are recurrently 
advanced with respect to the maxima given by the exact evolution. We also depict 
evolution under a non-hermitian effective Hamiltonian, and its effect is reflected in 
the fact that population maxima can be larger than one, as shown in the central max¬ 
imum depicted in red. Since we are looking at the low energy effective evolution, 
the fast oscillations do not appear in the evolution under the effective Hamiltonian. 

Analogously, in Fig. [3] one notices that the hermitian Hamiltonian to 10th order 
(obtained from computing the perturbative expansion of B to 10th order and substi¬ 
tuting in ( |2T| ) perfectly matches the evolution of the exact populations, when one 
discards the rapid fluctuations that populate temporarily the highest energy state. 
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Fig. 3 Evolution of the population of the ground and excited states, with initial state (1,0,0) 7 ^, 
under a) the exact Hamiltonian (continuous black line), b) h'^j- (dotted red line) and c) hy 0 ^ with 
V = 1 (dashed green line). The parameters are S = —0.01754, Q a = 0.44, = 0.34, as before. 


7 Large frequency expansions for periodic Hamiltonians 

Let us consider a time-dependent periodic Hamiltonian with period 2n/(0, so that 
H(t ) = H(t + 2n/a>). Let U{t) be the operator solution of id r U(t) = H(t)U(t), with 
U(0) = 1. Then, by Floquet’s theorem, the eigenvalues of U(2.71/co), when written 
as exp (— i27te/co), provide us with the so-called quasi-energies e. It is well-known 
that the quasi-energies can be identified with the eigenvalues of the corresponding 
Floquet’s Hamiltonian fl4l . In order to understand it, let us assume that the original 
periodic Hamiltonian can be expanded as 

H(t)=^H k e-‘ kcot (29) 

k 

with integer k and H k = // i-. Let us define the operators K and T acting on a 
different Hilbert space (the canonical example is L 2 (.S' 1 j), and such that a) K is 
hermitian and its spectrum consists of all integer numbers, i.e. if we denote the 
eigenbasis of K as |£), then K\k) = k\k)\ b) T is unitary; c) [K, T] = T. The Floquet’s 
Hamiltonian associated to the original one is hence 



( 30 ) 




Beyond Adiabatic Elimination: Effective Hamiltonians and singular perturbation 


13 


In order to avoid cluttering formulae unnecessarily we frequently omit the tensor 
product sign; by H^T k we actually denote Hk®T k , which acts on if 

is the original Hilbert space on which H{t) acts. Similarly, K stands for 1 0 K . 

We shall now apply the formalism of adiabatic expansions to this Floquet’s 
Hamiltonian to obtain approximate values for the quasi-energies. The relevant pro¬ 
jectors are P = 1 ® |0)(0| and Q = 1 — P. Therefore, 


PH f P = H 0 , (31) 

PH f Q = J]H k PT k Q = £ fl*® l°X-*l - (32) 

k t-^0 

QHf P = Yj H kQ Tkp ='E H k® \k) <0|, (33) 

k kyk 0 

QHfQ = -coQKQ + Ys H kQT k Q • (34) 

k 


The problem in the direct application of the formalism presented above lies in the 
computation of the inverse restricted to the Q subspace. However, in the limit of co 
much larger than H k in norm, the inverse may be approximated by 


Q(QHfQ)~ 1 Q 


v 


oo 1 

S 


Kg'Y i H t Q rk Q I K, 1 


Q ■ 


(35) 


where Kg 1 is the inverse of K restricted to the complement of |0). By using now 
this approximation and truncating, we obtain a first approximation of the effective 
hermitian Hamiltonian in the constant sector as 

H 0 --YjH_ k H k . (36) 

® k% k 

As an example admitting an exact solution, let us consider 

H(t) = g (a V®' + . (37) 

The effective Hamiltonian reads 

(38) 

co \ ® 2 / 


up to 0(g 6 /co 5 ). It should be stressed at this point that the crucial information car¬ 
ried by the effective Hamiltonian is the set of eigenvalues. As stated, this example 
is exactly solvable, but it is also an interesting toy model to test the power of the 
techniques shown in this paper. To solve it, it is sensible to consider the evolution 
of a time-dependent Hamiltonian, which is, in physics language, the interaction pic- 
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ture Hamiltonian obtained from Hs = (A + 0))(7 Z /2 + gO x with respect to a free 
part given by //f ree = coo z /2. For the sake of clarity, let us define precisely what the 
interaction picture is. Let Ugit) be the unitary solution of the initial value problem 
idfUs = HsUs with Us{ 0) = 1. Similarly, let t/f ree (f) be the unitary solution of the 
initial value problem id t Uf lee = //fr ee t^free with t/f ree (0) = 1. The interaction picture 
Hamiltonian with respect to the free part //f ree is defined as 

Hl(t) = [/ fr ee(0 f (H s - //free) t/f ree • 

Now, let us consider the initial value problem id t Uj(t) = Hj(t)Uj(t), with Uj( 0) = 
1 and Uj(t) unitary. Then one readily sees that the solution for this initial value 
problem can be written as 


u,(t) = Uf ree (tYu s (t). 


(39) 


Alternatively, if we are presented with a time dependent Hamiltonian, the corre¬ 
sponding evolution operator can be computed if the Hamiltonian is identified as the 
interaction picture Hamiltonian with respect to some free part. 

For the specific /A, and //f ree above we have 


// 7 (f) = e imaZ / 2 (^a z +ga x ^j e i(0,aZ / 2 

= — a z +g (cos (cot) a x — sin(cot)a y ) 
= ^(7 Z + g (o + e l(01 + cr~e~'® r ) . 


Obviously, the time dependent Hamiltonian H(t) of (37 i corresponds to this 
///(/) in the case A = 0. Particularising to the example the solution (39 1 we have 


.... ,i(ot 

///(/)= exp ( — o~ 


exp 


i((o + A)t 


o z — igta' 


(40) 


Hence, by computing Uj(2n/co), its eigenvalues read (2n/a )) x ^<b/2± \J{(0 -I- A) 2 /4 + g 2 
Thus, the quasi-energies, which are defined modulo co, can be fixed as 


=F 



( 41 ) 


If CO A,g, the quasi-energies may be expanded to obtain 



g 2 A g 2 (A 2 — g 2 ) 
CO 2 to 3 


matching (38 i in the case A = 0. 
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One can also apply Schrieffer-Wolff’s method in this case; no resummation is 
gained, however, since A ?\3 is computed exactly as 

^o = -|(ct + ®|- 1 )( 0 |-^®| 1 )( 0 |) , (42) 

giving = tjCH for the resonant case A = 0. The case A / 0 can also be solved 
exactly for resulting in 


^0= Z ^(a + ® | - 1)(0| -CT^ |1)(o1 ). 


(43) 


Therefore, in this case 4^0, the corresponding Hamiltonian is 




2^ 2 


A-2(0 


a . 


(44) 


To summarise, we have shown that both Bloch’s expansion and Schrieffer-Wolff’s 
provide us with large frequency expansions for the quasi-energies in periodic Hamil¬ 
tonians. As a further scope, it would be interesting to consider the application of 
these novel techniques to the quantum Rabi model beyond rotating-wave approxi¬ 
mation Gsmu. 


8 Conclusions 

The well established adiabatic elimination procedure in quantum optics has pro¬ 
duced some controversy in that scientific community on its meaning and on the fea¬ 
sibility of systematic improvements on that approximation. We have shown that, in 
fact, adiabatic elimination is the first term of a systematic expansion (be it perturba¬ 
tive or iterative) using Bloch’s equation from nuclear optics. We have further shown 
the existence of relevant solutions. We next connected this approach to Schrieffer- 
Wolff’s method, showing in which sense it can be said that Schrieffer-Wolff is a 
resummation. After the presented example, we have looked into the case of periodic 
Hamiltonians, employing the previous techniques on Floquet’s Hamiltonian, to ob¬ 
tain high frequency expansions for quasi-energies. We expect these novel methods 
to be useful in current models of interest as is the case of the quantum Rabi model. 
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